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Abstract 

We study the electric properties of single-stranded DNA molecules witli hairpin-like shapes in 
the presence of a magnetic flux. It is shown that the current amplitude can be modulated by 
the applied field. The details of the electric response strongly depend on the twist angles. The 
current exhibits periodicity for geometries where the flux through the plaquettes of the ladder can 
be cancelled pairwise [commensurate twist). Further twisting the geometry and changing its length 
causes complex aperiodic oscillations. We also study persistent currents: They reduce to simple 
harmonic oscillations if the system is commensurate, otherwise deviations occur due to the existence 
of closed paths leading to a washboard shape. 
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Introduction. — Recently, conduction properties of molecules have been intensely investi- 
gated. This has been fueled by a strong motivation to replace semiconductors with molecules 
that can save manufacturing expenditure due to their self-assembly properties. One of the 
significant steps towards molecule-based electronic devices has been made through DNA 
conduction measurements. Superconductivity was observed, but was supposed to be due to 
the proximity effect of the superconducting leads jlj. A random base sequenced DNA was 
shown to behave as an insulator P|, whereas the homogeneous poly(G)-poly(C) molecule 
have a large gap in their current-voltage (IV) characteristics Also, Ohmic behavior of 
DNA ropes was studied by Fink et al. 

On the other hand, in small systems where the phase coherence of electrons can be main- 
tained, quantum interferences play a key role in determining the characteristic properties. 
Indeed, a large number of devices that have been proposed, such as switches and transistors, 
are based on the wave nature of electrons Q|. One way to observe the interferences in a 
controllable way is to apply a magnetic field. It comes into play when the systems can 
accommodate non-simply-connected paths; the presence of the fiux tube gives a phase shift 
in the wave packet of the particle and changes the interference pattern j^. 

Despite abundant literature on DNA conduction, less attention has been paid to the 
existence of interference effects in this system. Molecules however have a huge variety in 
their structures allowing for different interference events. For instance, in a double stranded 
DNA, base pairs wind about the helical axis, which can be modelled as electrons travelling 
on a twisted ladder-like structure. When a magnetic fiux is present, there exist trajectories 
enclosing a finite fiux, which may affect the physical properties of the system. Furthermore, 
if one considers single stranded DNA or RNA, one may observe even more interesting ge- 
ometries with loops and bulges. The purpose of this paper is to elucidate the effects of 
interference on the electronic properties of these molecules. As a minimal model, we con- 
sider a hairpin-like structure in the presence of a magnetic fiux, as depicted in Fig. 1. The 
bases are coupled to each other via electron hopping and Coulomb repulsion. For simplicity, 
we consider them as identical except for the hopping constants. We evaluate the currents 
through the structures when they are coupled to electron reservoirs. The magnetic fiux 
applied to the system is shown to modify the current amplitude. It is also shown that the 
characteristics depend on the twist angles. For a geometry where the fiuxes through the 
plaquettes of the ladder are cancelled pairwise (we refer to it as commensurate ladder), the 
current oscillates as a function of the fiux through the loop in units of fiux quantum /, with 
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period / = 1, as in a single-loop Twisting the geometry or changing its length changes 
the periodic oscillation into complex aperiodic ones. We have also evaluated the persistent 
currents (PC) of the system which with a commensurate ladder behaves as that of a 
single loop. For an incommensurate ladder, non-vanishing contributions are built-up by the 
closed paths embracing the ladder, leading to washboard-shaped PC. 

System. — We assume that the system is described by a Hubbard-type model for spinless 
fermions jsj: 

H = ^[tijclcj + UijUiUj + h.c], (1) 

where rii = c\ci with Cj the operator annihilating the Fermionic particle at site and 
the sum is over the nearest neighboring pairs. This effective Hamiltonian was used to 

Q 

explore the conductance gap observed in the experiments of ref.j3| and the gapless states 
in the engineered DNA Partitioning the hairpin into C(hain), L(adder), and R(ing) 

sectors (see Fig. 1), we assign the hopping parameters along the contour as follows: 

h for 2-j = ±l,(z,j) e C 

t\\ exp(±27ri/) for i - j = ±1, (z, j) G R 
t|| exp(±27ri0j) for i — j = ±1, (hj) G L 




FIG. 1: Schematic picture of a hairpin-structured single-stranded DNA where bases are paired in 
the L (adder) sector, unpaired in the R(ing) and in the side C(hains). The right figures are shown 
to illustrate effective-flux configurations in the ladder with the twisted rung-bonds. 
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In addition, the ring hoppings in the ladder are characterized by 

ti<j ^t±ex.p{2nif) (2) 

and ti>j = t*^j with i and j satisfying i + j = Nc + Nl + Nr + 1 with Nc,l,r the number 
of sites in the chain, the ladder, and the ring, respectively. 

Here the phases accompanying the hopping integral are due to the presence of the thread- 
ing flux / and are defined as Uj — i||e*'^*^ with Aij — (27r/$o) Ji ' d£. For our system, 
the phase factor in the ring takes a simple form as Aij — 27r$/7VR$o = 27r//A^R with $ 
measuring the total flux through the ring in units of the flux quantum $o- The magnetic 
flux penetrating the unit cell in the ladder introduces the phase factor 9i given by 

Oi = S.Bay % (3) 

where a is the effective path length between bases, is the factor accounting for the area 
variation along the sites due to the twist of base: 



and 27r/T is the twist angle that can be varied from sample to sample, for example, by 
changing the amount of added gyrase during preparation. Relating the magnetic fluxes 
through the loop and the unit cell of the (flattened) ladder through the area ratio 7 = 
A7t/{Nr + 2)2, we deflne 9. = ^Sif. 

Let us consider the case of strong Coulomb repulsion where the ground state of the system 
at half-fllling contains charge density waves {ni) — (l-|-Acos(7ri-|-'0))/2. At low temperature 
(U/kBT ^ 1), one can safely consider the mean-fleld version of the Hamiltonian 

y^MF ^ "^[f.-clcj + h.c] +UY1 ^^^^^^ + 

(hi) i 

where the renormalized strength of the repulsion is given by C/ — U/S.. Therein, we have 
made the asumption that the repulsion between bases in pair in the ladder is negligible and 
ip is constant over the system. If this is not the case, an intricate point arises; specifically, 
when Nr is odd, the density waves must be distorted and have abrupt spatial changes in 
■0, giving rise to a kink (or a soliton). The existence of kinks and their dynamics in an 
odd-numbered ring can be another problem of interest. However, for the sake of simplicity, 
we leave it for future studies. 
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Transmission and currents — Let us now study the physical properties of the system by 
evaluating its electric conductance. To this end, we consider electrodes coupled to the system 
as in Fig. 1, introducing the self-energy correction Sjg^^^ on the edges of the chain as 

= ^lead(^i^i + "^iv^^)- (^) 

There is in fact some energy dependence in Sjg^^^^ but considering the bulky electrodes mostly 
used in experiments, we regard it as energy independent and identical for both electrodes. 
Evaluation of the Green function of the system, G{E) = {E — — He — /0+)^^, leads to 
the transmission coefficient t from the left to the right lead as 

t{E) = TrTeG{E)T, (7) 



where the coupling matrices are given by r^(r) = lm{TiYQg^(^)Sg(^r),i{N)- The corresponc 
currents driven by the finite voltage bias are obtained by the Landau-Biittiker formula 
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dE[fiE-^^,)-fiE-^^r)]\t{E)\', (8) 
-CXD 

where f{E) is the Fermi-Dirac distribution function, H£ and /i^ are the chemical potential 
of the left and right leads respectively, the difference of which is controlled by the applied 
bias voltage V as fie — fir = eV. 

Figure 2 displays transmissions and currents for the two configurations (see the upper 
panel), the ladder and the hairpin with their ends coupled to external leads. The trans- 
mission windows for the hairpin gets narrower than that for the ladder, due to the smaller 
number of bonds. For the ladder, when magnetic fields are absent, the IV curve agrees with 
the IV characteristics for the homogeneous poly(G)-poly(C) molecules with a large gap j^. 
When a magnetic field with a strength / = 0.5 is applied, the current is suppressed. This is 
also the case for the hairpin, with a different degree of suppression. 

To see how the magnetic field alters the current amplitude, let us fix the voltage bias, for 
example, V = 0.3 eV and V = 0.4 eV for the hairpin configuration of Fig. 2. Figure 3 (a) 
demonstrates how the corresponding currents oscillate with varying fiux. Here we present 
three curves for a given voltage, each one for a different twist angle. The interesting feature 
is that for T = 10 (solid lines) the oscillation is periodic with period / = 1. On the 
other hand, oscillations are not complete even until / = 2 for the case T = 10.5 (dotted 
lines) and T = 11 (dot-dashed lines). One can also notice that for T = 10, the current 
dips are located at the multiples of / = 0.5. Taking other twist angles moves the dips. 



This twist-angle dependence can be understood as follows: Consider the area factor Si. 
The flux through a closed path cj in the ladder is given by Xliecj =~ ^^'^^ ■ ^or the 
hairpin configuration (Nq = 4:,Nl = 16, T = 10), every path Ck has its respective counter 
path Cj in such a way that pairwise cancellation of the flux occurs since 
Consequently, the flux through the ladder gives null effects. In that case, in terms of the 
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FIG. 2: The transmission vs energy (lower panel) and IV curves (upper panel) for the ladder and 
the hairpin (see the schematics in the upper panel). We take 30 base pairs for the ladder, and 
Nc = 4, Ni = 14, Nr = 8 for the hairpin, and T = 10. The band parameters used throughout this 
paper are = 0.2 eV, t± = 0.1 eV, U = 0.2 eV . The dotted curves are for / = and the solid 
curves for / = 0.5. 
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FIG. 3: The current vs magnetic flux at a given voltage bias for hairpins in different configurations 
with Nc = 4: (a) Nl = 14 and Nr = 8; (b) Nl = 30 and Nr = 14; (c) Nl = 30 and Nr = 26 with 
T = 10 (solid lines), T = 10.5 ( dotted lines), and T = 10.5 (dot-dashed lines). 



magnetic flux, the geometry is topologically equivalent to that of a single loop connected to 
the electron reservoirs (see Ref. |2||): The transmission of this configuration has been shown 
to oscillate with period / = 1. When twisting the structure slightly out of the commensurate 
configuration (for example, T = 10.5), such a perfect cancellation does not occur, leading 
to complex aperiodic oscillations (within the flux interval) . One can also expect that adding 
more plaquettes must yield similar features. In Fig. 3 (b) and (c), doubling the size of the 
ladder, we plot the currents at = 0.4 eV. Neither of the twist angles leads to periodicity. 
Also, the current enhancement by flnite flux < / < 0.5 can be noticed. Let us here give 
an order of magnitude of the strength of the magnetic field for the effects to be measurable. 
For small loops with radius 50 (corresponding approximately to 8 bases), in order to achieve 
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FIG. 4: The persistent currents vs flux: The upper panel shows the PC of the hairpins {Nc = 
4:,Nl = 14) with different number of bases in the loop, Nr = 6 (solid line), Nr = 8 (dashed line), 
Nr = 10 ( dot-dashed line), and Nji = 12 (dotted). While the oscillation is periodic in /, the 
amplitude decreases as the loop size increases. In the inset, we show that the curves coalesce when 
using a scaling IpC^'^^^^. The lower panel displays simple harmonic oscillations of PC for Nc = 4, 
Nr = 6 and = 14 (solid line), obtained by changing the ladder size = 18,22,26 ( dotted, 
dashed, dot-dashed, respectively). The inset shows the fits according to Eq. (|T2|l for Nl = 2 and 
iVL = 4. 

/ = 0.5, one needs fields in the range of thousands of Tesla. Considering that at present 
the maximum field strength is about 100 Tesla jisl, it would be easier to realize a larger 
loop. In Fig. 3 (c), we take 24 bases in the loop, making the area about ten times bigger 
than that in Fig. 3(a). Even if the field strength can be in the realizable range, the effects 
remain still significant and are not washed out by the enlarged loop-size. 



8 



Persistent currents — Now we evaluate the persistent currents when the system is isolated. 
This allows to assess how much the magnetic fields alter the energy levels (whereas the 
magnetic-flux effect demonstrated above mixes the influence on wavefunctions and on the 
energy spectrum). This becomes clear when writing the Green function in an expansion of 
eigenfunctions: 

leading to the transmission coefficient given in Eq. Here the magnetic fiux affects not 
only the energy spectrum En but also the wavefunction ip{x). Let us take a simple example 
for understanding persistent currents. Consider a ring of radius R threaded by a fiux </> 
(in units of the fiux quantum $o)- The angular momentum of free electrons traveling on 
the ring is shifted by the magnetic fiux as {n\L\n) = h{n + 0) with (n|^|n) denoting the 
expectation value of an operator A in the nth eigenstate. The energy eigenvalues are then 
given by £n = h^iji + (pY /2mR^. At zero-temperature, the evaluation of currents in the 
ground state is straightforward 

n n 

where v is the velocity operator, and the sum is over all the occupied levels, n = 
0, ±1, . . . , ±(A^e/2 — 1), —Ne/2 for an even number of electrons and n = 0, ±1, . . . , ±(A^e — 
l)/2 for odd Ng . It is clear that the ground state carries currents so that if dissipation is 
absent in the system, the currents keep fiowing, and are thus named persistent currents. 
The currents oscillate with the magnetic fiux with a period of = 1, and their ampli- 
tudes decrease as 1/R. When band gaps are present in the system (due for example to 



Coulomb repulsion or coupling to 
exponentially with increasing 



attice distorsion), the persistent currents are suppressed 
In this case, one has 



/p = (-l)^^/oe-"^sin(27r/) (11) 



with a = for vanishing gap. Here the sign of the current depends on the parity of 

Our system (a loop coupled to a ladder) has a number of closed paths or rings, with 
different areas and sizes. There are contributions of persistent currents from all the rings, 
and therefore, Ip = where /j denotes the current fiowing in the i-th ring of size Ri, 

enclosing the fiux /j. We first evaluate the persistent currents in the hairpin with Nc = 4: 
and Nl = 14 where the magnetic fiuxes through the ladder plaquette effectively cancel one 
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another (see the upper panel in Fig. 4). The currents behave as those of a single loop; note 
that the sign of the current changes with the parity of Ne (we consider half-filled systems) . 
Furthermore, increasing the number of bases decreases the current amplitude. Following the 
finite size scaling in Eq. (fTTll . all the curves are shown to be merged into a single one (see 
the inset). Adding more plaquettes causes a deviation from a single-component sinusoidal- 
curve (see the upper panel in Fig. 4). Since larger rings die out exponentially , taking paths 
enclosing up to two unit-cells, one obtains the approximate form of the currents as 

Ip ^ Ir sin(27r/) + ^ sin(27r0i) + ^ Id,i sin(27r0i), (12) 

i i 

where 6i is the fiux through the iih plaquette of the ladder, and (pi is the sum of two 
neighboring 6iS. For hairpins with a few plaquette in the ladder, Eq. fT^ perfectly fits the 
numerical evaluations (obviously, for a hairpin with a single plaquette, the above expression 
is exact). 

Summary and Remarks — We have demonstrated how a magnetic flux could influence 
the electronic properties of molecules with a twisted hairpin-like shape. The current has 
been shown to oscillate with flux changes. Geometry factors such as twist angles and the 
number of rungs of the ladder play a crucial role in determining oscillation patterns. We 
point out that there exists geometries where the flux through the ladder vanishes, so that 
the sole cause for the oscillation is the flux through the loop. For that case, the current of 
the ladder ressembles that of a ring. Twisting the geometry brings about complex aperiodic 
oscillations. These oscillations are also reflected on persistent currents. The non-vanishing 
contributions of the closed paths embracing the ladder leads to washboard-shaped persistent 
currents. 

We should note that the above mentioned results were predicted under the assumption 
that electrons preserve their initial phase in their propagation. Therefore, maintaining phase 
coherence is crucial in experimental observations of these phenomena. The phase coherence 
is strongly affected for example by temperature (or generally speaking, coupling to the 
environmental degrees of freedom). For DNA molecules for which conduction properties 
are measured at room-temperature, phase coherence can still be traced. However, finite- 
temperature effects cause fiuctuations in the structure. Even if the effect is small, the 
plaquette areas as well as the phases due to the fiux become randomized. This addresses 
the issue of (weak) localization that will be considered in a future work. 
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